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The interactions of liquid metals and alloys with the environment mostly depends on the 
thermodynamic properties of the liquid surface. In fact, the surface tension is strongly 
influenced by the presence in the surrounding atmosphere of reactive gases through 
solution, adsorption mechanisms and/or surface reactions. In particular, oxygen, which 
shows a high surface activity towards a large number of metallic systems, is the most 
important contaminant of liquid metals and alloys. 

Theoretical approaches for estimating the oxygen mass transfer at the liquid vapour 
interface under inert atmosphere and vacuum have been developed already in orderto relate 
the observed physical properties to the real surface composition data. 

In the present work a model of the interfacial transport of a liquid metal oxygen system 
under Knudsen conditions that foresees the temporal evolution of the interfacial 
composition is presented. The diffusion characteristic times for reaching steady-state 
conditions are evaluated in order to define two system "sizes" depending on the different 
oxygen transport mechanisms in the liquid phase. 

An experimental study of the interface evolution is at present under way and preliminary 
results show a satisfactory agreement with theoretical studies. 


Nomenclature 

c = molar concentration (mol m -3 ) 

D 0 = diffusion coefficient of atomic oxygen in the 
metal (m 2 s - 

/ = function defined in Equation 41 
k = partial transport coefficient for the mass transfer 
(mol sg -1 m -1 ) 

K = transport coefficient (m s -3 ) 

M = molecular weight (g mol -1 ) 

N = molar flux between the gas and the condensed 
phase (mol m -2 s - J ) 

Pj = partial pressure of the j th oxide (Pa) 

P M = vapour pressure of the pure metal (Pa) 

P G = partial pressure of molecular oxygen (Pa) 

P lot = total pressure (Pa) 
r = radial co-ordinate (m) 

R = drop radius (m) 

M = gas constant (J mol -1 K - x ) 
s = liquid saturation degree 
T = temperature of the system (K) 

T b = boiling temperature (K) 

T m = melting temperature (K) 

T w = temperature of the container wall (K) 
f = time (s) 

x 0 = molar fraction of atomic oxygen in the liquid 
phase 


y = molar fraction in the gas phase 
oc = condensation coefficient 
r = Gibbs adsorption (molm -2 ) 
v = stoichiometric coefficient 


Subscript and superscript 

b = bulk 
D = diffusion 
1 = liquid 

g = gas 

0 2 = molecular oxygen 
M = metal 
j =f h oxide 
I = interface 
0 = initial state 
s = saturation state 
st = steady state 

1. Introduction 

The studies on the interplay of all possible mechan¬ 
isms of mass transport at the interface of liquid 
metal-gas systems are of a topical interest, from a sci¬ 
entific and a technological point of view. In particular, 
for the high temperature processes, the definition of 
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the factors controlling oxygen mass transport through 
the liquid metal interface as a function of the bound¬ 
ary conditions assume an important role also in the 
definition of the efficiency of technological processes. 
In fact, oxygen shows a strong surface activity towards 
a large number of metallic systems [1], and its pres¬ 
ence, in practice, can never be avoided. For this rea¬ 
son, investigations on high-temperature capillarity are 
basic to evaluate the tensio-active effects of oxygen on 
the liquid metal surfaces through surface tension and 
contact angle measurements. 

Because the high temperature capillarity processes 
are generally performed under controlled atmosphere 
or under a vacuum, it seems extremely useful to pro¬ 
vide theoretical models, in addition to the experimen¬ 
tal observations, which enable an accurate estimation 
to be made of both the degree of contamination of the 
surface, and the mechanism of the oxygen mass trans¬ 
fer at the liquid-gas interface, in order to relate the 
observed physical properties to the real surface com¬ 
position data. 

To attempt this goal, we have previously developed 
two theoretical descriptions [2, 3], on the basis of 
kinetic and transport theories, relating the mass ex¬ 
change between liquid metals and the surrounding 
atmosphere at different oxygen partial pressures: the 
first one under inert atmosphere (He fluxing with an 
inlet pressure of % 10 5 Pa); the second one under 
a vacuum (total pressure lower than 1 Pa in the 
Knudsen regime). 

In the first model, where diffusion is the main mass- 
transport mechanism in the gas phase, reactions 
between metal vapours and molecular oxygen are 
allowed to occur, establishing specific reaction re¬ 
gimes in which the oxygen can or cannot contaminate 
the liquid metal surface. 

In the second one, a total pressure less than 1 Pa is 
considered, assuring the Hertz-Knudsen conditions of 
thermomolecular regime for the mass transport phe¬ 
nomena in the gas phase. In fact, under a vacuum, 
transport processes in the gas phase have no control 
over the total interphase mass exchange, because of 
the large mean free path of the molecules. 

Moreover, when dealing with metals which form 
volatile oxides, condensation of the oxide vapours on 
the reactor walls must be taken into account; the 
removal of oxide vapours from above the metal sur¬ 
face can significantly displace the oxidation equilib¬ 
rium, enhancing the rate of oxide vaporization; this 
effect can be exploited in some cases to deoxidize 
and/or purify the liquid metal. 

Owing to the double contribution of molecular oxy¬ 
gen and volatile oxides to the oxygen flux to and from 
the surface, steady-state conditions can be established, 
in which no variations of composition of the two 
phases occur over time. The problem can be ap¬ 
proached by evaluating the total oxygen and metal 
evaporation rates as functions of the overall thermo¬ 
dynamic driving forces. 

Curves separating different “oxidation regimes”, 
were drawn, which represent an accurate control 
tool of the oxygen contamination degree of liquid 
metals under a vacuum; they are useful in capillarity 


measurements, in order to relate the observed surface 
properties to the actual liquid composition. 

In this paper we present a theoretical model, based 
on the assumption of the model described above and 
upgraded, concerning the oxygen transfer through the 
liquid-vapour interface under Hertz-Knudsen condi¬ 
tions, allowing the temporal evolution of the inter¬ 
facial composition to be considered. The model allows 
the diffusion characteristic times for reaching steady- 
state conditions to be evaluated. 


2. The model 

Let us consider a liquid metal drop placed at the 
centre of a close chamber as an exemplification of the 
real experimental set-up. The following conditions 
apply: 

1. The temperature near the drop is maintained at 
a uniform and constant value T, which is higher than 
that on the chamber walls, T w . 

2. Total pressure P tot < 1 Pa is maintained inside 
the chamber: thus, Hertz-Knudsen conditions are as¬ 
sured for the mass transport in the gas phase [4]. 

3. The activity coefficient of the liquid metal is 
assumed to be unity. 

4. Small amounts of dissolved oxygen are assumed 
not to affect significantly the vapour pressure of the 
pure metal, P M , which is assumed to depend only on 
the temperature T. 

5. The partial pressure of oxygen in the residual 
atmosphere, P 0 ,. g is considered to be constant. 

6. A molar fraction x 0 < x 0 , s of atomic oxygen is 
dissolved into the metal, where x 0 . s indicates the limit 
of solubility [5], Above this limit, the formation reac¬ 
tion of the first stable oxide in the condensed phase: 
M (cond) + vOf s) = MO ( 2 “ nd) , is shifted to the right. 

The term "first oxide” denotes the oxide with the 
minimum oxygen content which is thermodynami¬ 
cally stable under the imposed conditions. 

When x Q = x 0 . s , that is when the metal and the 
oxide exist at the same time as different condensed 
phases, the atmospheric composition is fixed at con¬ 
stant T for a system at equilibrium, and it can be 
calculated from thermodynamic data. In particular, 
the vapour pressures P j s of the species containing 
oxygen (suboxides) resulting from thermal decomposi¬ 
tion of the oxide depend exclusively on T [3]. 

The oxygen diffusion inside the liquid metal drop 
and the interfacial interactions can be evaluated as¬ 
suming the liquid metal drop as a rigid sphere. Under 
this assumption and under the conditions described 
above, the local oxygen mass balance in the drop can 
be expressed in spherical co-ordinates as 


0X ° n 1 s f 2 Sx o\ 

iV Do ?s( r '^j 


0 < r < R(t) < R 0 


( 1 ) 

where r is the radial co-ordinate which is zero at the 
centre of the drop, D 0 is the diffusion coefficient of 
atomic oxygen in the metal, R(t) and R 0 are the sphere 
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radius at time t and at initial time t = 0, respectively. 
The possibility of a drop shrinking, due to the metal 
evaporation, is explicitly taken into account. 

As initial conditions (t = 0 and R = R 0 ), we assume 
a uniform distribution of oxygen inside the drop 

x 0 = Xo for 0 < r < R (2) 

The assumed boundary conditions are: the symmetry 
condition at the drop centre 

f > 0 r = 0 ^ = 0 (3) 

or 

the oxygen balance at the moving interface 
f > 0 r = R(t) 

_ Cl ^°~dr — 2N ^ = ^ ClX0,1 — c s^ko 2 ,i) (4) 


4), the gradient of x n at the interface can be defined 


0x o 

“07 


— 21Y 0 ,i + 



2ko 2 , 


x (2iVo 2l i + 



(DqCi) 1 


(7) 


As Equation 7 shows, the gradient of oxygen molar 
fraction depends on the oxygen flux at the gas-liquid 
interface, on the concentration discontinuity between 
the interfacial gas and the interfacial liquid and on the 
interface displacement which is correlated to the oxy¬ 
gen and metal fluxes. 

If the condition c g 2y 0 , , I «c 1 x 0 ,i is assumed, which is 
a more restrictive condition than CgCccj but generally 
verified because of the rarefaction of the gas at the 
total pressure considered by the initial condition, the 
expression of the gradient of oxygen molar fraction at 
the interface (Equation 7) yields 


In Equation 4, c x and c g are the molar concentrations 
of the liquid and of the gas, respectively, N 0l l is the 
flux of molecular oxygen at the interface, y 0 ,,i and 
x 0 ,i are the molar fractions of molecular oxygen in the 
gas and of atomic oxygen in the metal res¬ 
pectively, near the interface. Equation 4 shows the 
oxygen flux arriving from the bulk to the interface 
( — CiD o (0x o /0r)) and the oxygen flux going through 
the interface towards the gas phase {N 0 ,,i), to be equal 
to the oxygen accumulation at the moving boundary. 

The variation with time of the surface adsorption 
(dT/df) is not considered in Equation 4; in fact, with 
reference to the known values of T for oxygen in liquid 
metals [1] and to the drop dimension considered 
(R = 5 x 10“ 3 m), the ratio between the oxygen mole 
number in the liquid bulk and those at the interface 
results much greater than one: i.e. (_RciX 0 /3T)» 1 (for 
instance of the order of 10 3 ). 

In addition to the boundary conditions, the total 
mass balance of the system has to be considered in 
order to deduce the interface evolution as a function of 
the net evaporation fluxes, both the oxygen flux at the 
interface iV 0 ,,i and the metal flux N MA coming from the 
interface 


0X O 

dr 


[(1 — x 0 ,i)2N 0j j — x 0 ,i N m ^(DqCi) 1 


( 8 ) 


As a further approximation, when x 0 .i«l is as¬ 
sumed, which is reasonable in many cases, the gradi¬ 
ent at the interface simply becomes 


0x o _ 2N 0i j 
dr , D 0 Ci 


(9) 


Under this approximation, the oxygen concentra¬ 
tion gradient at the interface depends on the evapor¬ 
ation flux and on the oxygen diffusivity ratio, but is 
independent of the evolution of the interface. 

Taking into account Equation 9, the oxygen bal¬ 
ance expressed by Equation 4 becomes 


d.R 

(cjxo.i - c g 2y 0 ,.i) — = 0 


which, on the basis of the assumption made above, can 
be interpreted also as 


d R 
d t 


x 0 


( 10 ) 


— (2N 0 ,.i + N M:l ) — [cj 


c g(2)’o,,i 


+ Vm, i)] 


d.R 

dr 


which confirms that the variation in size of the drop is 
negligible. 


(5) 

where Vm.i is the molar fraction of the metal in the gas 
phase at the interface. 

Because the concentration (density) of the gaseous 
phase is negligible compared with that of the liquid 
phase (i.e. CgCccfl, the total mass balance (Equation 5), 
can be written in a more simple form, deducing the 
interface evolution which is directly dependent on the 
net evaporation fluxes 

dH__2Ao,,i + N M j 

dr Cj 


3. Evaluation of steady-state conditions 

The steady-state conditions of the system corresponds 
to a uniform and constant drop composition 

x 0 = x OI and = 0 for 0 < r < R (11) 


In particular, the liquid-vapour interface is con¬ 
sidered to be in steady state when the oxygen concen¬ 
tration gradient in the liquid phase is zero, i.e. 


dx„ 

“07 


= 0 


( 12 ) 


Substituting the total mass balance (Equation 6) in 
the expression of the total oxygen balance (Equation 


As the condition to reach the steady state at the 
liquid-vapour interface (Equation 12) is implicit in the 
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condition of the steady-state of the system (Equation 
if), it can be stated that the steady-state of 
liquid-vapour interface is a necessary condition to 
reach the steady state of the drop. Obviously, such 
a condition can be verified when all fluxes at the 
liquid-vapour interface are zero 

No 2 ,i = N M ,t = 0 (13) 

Equation 13 defines the thermodynamic equilibrium 
at the liquid-vapour interface which is consistent with 
a null gradient and a steady interface, and results in 
the kinematics steady-state condition dR/dr = 0. 
Elowever, the steady state at the interface can also be 
achieved under less restrictive conditions than that 
imposed by Equation 13; in fact Equation 12 is 
verified (see Equation 8) also when a sort of flux con¬ 
gruence is achieved 


Equations 2, 3 and 9, obtaining 

dx 0 _ 6 A^o,,i 

dt Rc l 


( 19 ) 


Then, with reference to the mean value of the flux 


N n 


(N 


o 

o 2 ,i 


+ Nt.i) 
2 



( 20 ) 


(where N o,.i and N g 2iI , are the interface oxygen fluxes 
at the initial state and at the steady state, respectively), 
and to the initial (x 0 = xg) and final (x G = xg) com¬ 
position, the approximate time integration of Equa¬ 
tion 19 allows the characteristic time to be evaluated 




1 Rc i ,_.o 


3 N° 


(xg - xg) 


o,,i 


( 21 ) 


-fVo,,i _ x 0 .i 
Nm,i 1 — x 0 ,i 


(14) 


In fact, the presence of non-zero fluxes does not neces¬ 
sarily modify the composition of the liquid phase; on 
the other hand, when a net flux of matter crosses the 
liquid-vapour interface, under the steady condition 
(Equation 14), a time variation of the drop dimensions 
occurs 


d R 2N 0l ,i 

d t x 0 ,iCi 

(15) 

Moreover, under the frequently verified condition 

for which the oxygen flux at the interface is 
the metal evaporation flux 

less than 

2!V 0 ,,i«1Vm,i 

(16) 

Equation 14 can be simplified to 


Af 0 ,,i = 0 and A\ri ^ 0 

(17) 


This approximation is consistent with Equations 9 
and 10. In this case the velocity of the interface can be 
evaluated from Equation 15 while the study of 
transitions to the steady-state condition can be carried 
out under the approximation of constant dimensions 
as expressed by Equation 10. 

In the following, we refer to this latter approach, 
which seems to be sufficient to approximate the ex¬ 
perimental conditions in the short times. 


By comparing the characteristic times f D , as defined 
in Equation 18, and f st , two subclasses of systems can 
be distinguished: “small” and “large” drops, depend¬ 
ing on the different mechanism of the oxygen trans¬ 
port in the liquid phase. 


4.1. "Small" drops 

We define the drop as “small”, when the diffusion time 
within the drop is very small (related to the system 
parameters), so that the drop composition can be 
taken as uniform (x 0 = x OJ ; 0 <r < R',t > 0) and the 
condition f D «f st is verified. From a geometrical point 
of view, taking into account Equations 18 and 21, this 
condition is equivalent to 

R« 3 ^(x g-xg) (22) 

o 2 ,i 

As the “small” drop, defined in Equation 22, shows 
uniform composition, the characteristic time of the 
bulk, f b , and that of the interface, fi, are equal to f st , 
then 


t\ — t b — t sl y>ty, 

Thus, for “small” drops, the characteristic time is 
calculated according to Equation 21 and the steady- 
state condition x 0 = x 0 ,i = xg = const, is achieved at 
the same time both at the interface and in the bulk. 


4. Characteristic times 

In order to examine the physical meaning of the 
steady-state condition and to evaluate the time neces¬ 
sary to reach it, some characteristic times have to be 
evaluated. The diffusion time within the drop is de¬ 
fined [6] as 


The time f st necessary for a homogenous drop to reach 
the steady-state composition xg can be estimated first 
by a space integration of Equation 1 together with 


4.2. "Large" drops 

We define the drop as “large” when r D »t st and Equa¬ 
tion 22 becomes 

R » (xg- xg) (22a) 

•W o,,i 

In this case the diffusion within the liquid metal drop 
is the controlling process. 

The time f b is of the same order of r D , while should 
be much shorter and could be evaluated equalling the 
diffusional penetration distances. 

The same result can be obtained by following a for¬ 
mal approach, allowing the interfacial composition to 
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be calculated from the convolution integral 


* 0,1 


Ci(nD 0 ) 112 


(r -1') 1/2 


df (23) 


Equation 23 is correct for planar geometry, but it is 
fully acceptable also for spherical geometry if the 
diffusion penetration distance is much smaller than 
the drop radius. 

Substituting the mean value of flux (Equation 20) in 
Equation 23, an estimation of the characteristic time 
of the interface q, is obtained 


t, = n 


D 0 c, (xo — xg ) 2 

(2A 


(24) 


A correlation between the characteristic times can 
be obtained by using Equations 18 and 21 and assum¬ 
ing the approximation: n 1/2 x 2 


q = 


1 St 


(25) 


So, for “large” drops: q = t 2 t /t D «t st «t b ~ q> In this 
case the interface reaches the steady-state composition 
in a characteristic time q defined by Equation 24 or 25, 
which is much less than f b . This last time is also much 
longer than the time q, associated with a homogene¬ 
ous drop, and it is practically controlled by the diffu¬ 
sion time f D . 


5. Temporal evolution of the interface 

The trend of the oxygen interfacial composition ver¬ 
sus time, the initial conditions defined in Equation 
2 being valid, can be described by the differential 
Equation 19 for “small” drops, and by the integral 
Equation 23 for “large” drops. The inequalities (Equa¬ 
tions 22 and 22a) discriminate between “small” and 
“large” drops. 

In particular, as a simple case, the example of an 
overall oxygen flux linearly dependent on the driving 
force is reported. By assuming for the overall atomic 
oxygen flux 

2!Vo 2i i = Kci(x 0 i — Xq) (26) 

(where K is the transport coefficient for the oxygen), 
the characteristic times of the interface for the two 
subclasses considered, on the basis of Equations 21 
and 24 are, respectively 

2 R 

q = - — for “small” drops (27) 

3 K 

and 

4 D 0 

q = —y f° r “large” drops (28) 

K 


On the other hand, for “large” drops, also Equation 
23 admits the analytical solution 


Xq.i - Xq 
Xq - Xq 


= e 4t/t| erfc(4t/q) 1/2 


(30) 


In particular, for t»q Equation 30 degenerates to 


Xq.i - Xq ^ / q y /2 ^ J_ / DoV 12 
Xo — Xq \4jrfy K\ntJ 


(31) 


which is an expression of an approximate trend of the 
dependence on time of the oxygen interface composi¬ 
tion in the liquid phase for “large” drops. In this case 
also, the expression of the flux at the interface degener¬ 
ates in the well known trend at constant composition 
*o,i 

/ n \ 1/2 

2!Vi * (xg - xg) (32) 


However, in more general cases, where the inter¬ 
facial flux Ai 0 ,.i is non-linear, the characteristic times 
also depend on the concentrations. 

As already described in the previous model [3], 
concerning the oxygen mass transfer at liquid- 
metal-vapour interface under a low total pressure, 
with its assumptions remaining valid, the total oxygen 
and metal fluxes at the interface are written as 


No,,I — k o 2 {s 2 Po 2 ,s — Po 2 ,g) 


+ XjVjkj(s?'Pj' t - P Jig ) 

(33) 

Nm.i = k M (P Ms — PM.g) 


+ Xjkj(s?‘Pj, s - Pj g ) 

(34) 


where k Qi , k M , and kj, are the coefficients of the Hertz 
Langmuir-Knudsen [4] equation for the evaporation 
from a surface under a vacuum. For an evaporating 
element or a compound, such a coefficient is given by 

( 1 y/ 2 

k = a (-= ) (35) 

\2n.MT M ) y 

where a is the condensation coefficient of the evapor¬ 
ating species and M is the molecular weight. 

Equations 33 and 34 describe the system in terms 
of coefficients and driving forces related to the gas¬ 
eous phase only. In the expression of the driving 
forces appears the relative saturation of the interfacial 
liquid. 

The steady-state condition for the interface (Equa¬ 
tion 14 or 17) combined with Equations 33 and 34 
leads to individuate the value of the steady-state liquid 
saturation degree s st , as a function of the gaseous 
phase composition (P 0 ,, g - Pj . g , F Mg ), of the values at 
saturation (P 0 ,, s , Pj.s, F MiS ) and of the liquid phase 
(x 0 , s ), which in turn depend on temperature. In the 
expression of s st also the transport coefficient ratio 


For “small” drops, the analytical solution of 
Equation 19 yields an exponential dependence on 
time of the oxygen interface composition in the liquid 
phase 


Xq.i — xg 2t/t , 
xg - xg 


(29) 


k j _ 1 

k 0l a 0i (m + Vj) U2 


(36) 


appears, where m is the molecular weight ratio = 
Mm/Mq,. By taking into account Equations 17, 33 and 
36, the “steady-state saturation degree” s st = xg/x 0>s 
can be derived by the following equation, which is 
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(2) if s 0 « s st 


solvable by an iterative procedure 

1/2 

(37) 

As expected, when the saturation values of oxygen 
pressure, P 0 ,. g = Po 2 , s (T) and the vapour pressure of 
species containing oxygen, P js = Pj, s (T ) are imposed 
(oxygen saturation, with the system in monovariant 
equilibrium), s st becomes equal to unity. 

When, as in this case, the steady-state condition can 
be approximated with Equation 17 (1 Vo 2 ,i = 0), the 
oxygen flux at the interface (Equation 33) in the steady 
state can be written in the form 

A1o 2 ,i = ^o z ,\Po 2 ,s( s i ~ S S) 

+ XjVjkjMs?’ - s;,') (38) 

and, by the same way, the initial flux N o„i, which is 
non-linear with respect to the term (so — sj), results 

Aio 2 ,I = ko,,lPo 2 ,s(So — S st) 

+ ZjVjkjjPjJsp - .s;, v ) (39) 


s s , = 


• O;. g A ■ j. g 

+ ZMjPj.sS^' " I 


Then, the characteristic time f st , as defined in Equa¬ 
tion 21, depends on the initial composition of the 
liquid and, on the basis of Equation 39, it can be 
written as 


fst 


RciXq.. 


1 


3^o 2 Pq 2 ,s /( s o> s st ) 


where / is 
/='" 2 " 2 


si-sl+XjVjh-^isp-s?) 

M) 2 ,s 


(40) 

(so s st ) 

(41) 


6. Applications 

For the sake of simplicity, the factor / can be approx¬ 
imated under different conditions as follows: 


(1) if s 0 »s st 


So 


Pi 

Po 2 


i +¥A^^ (vrl) 


(42) 


/= 2s 0 


l+S^^-^-sg^ 11 

r 0 2 ,s 


(43) 


in particular, Equation 42 can be usefully applied to 
calculate t st when the initial condition of the system is 
the saturation (s 0 = 1)- In this case, the condition 
s 0 »s st is largely verified and the factor/can be further 
approximated in 


/ 


1 


+ ZjVjh 



(44) 


Under these conditions the characteristic time f st can 
be calculated by the relationship obtained from 
Equations 40 and 44 


fst 


Rc\X q,s \ 

3ko 2 Po 2 ,s/ 


1 + XjVjfy 



(45) 


The characteristic time t st for some metallic systems 
of technological interest has been calculated by using 
Equation 45, considering, as a first approximation, the 
ratio of the condensation coefficients oq and a Q , to be 
equal to unity if the efficiency of the evaporating 
process is assumed to have about the same value for 
both species. 

The thermodynamic quantities have been cal¬ 
culated as a function of temperature starting from the 
melting temperature up to the boiling temperature at 
the fixed total pressure of 10“ 3 Pa. The partial pres¬ 
sures P 0 , s and Pj' S have been calculated by evaluating 
the equilibrium constants of the reaction of formation 
of the first stable oxide and the volatile oxides, respec¬ 
tively [7], on the basis of equilibrium constant of the 
oxidation reaction. 

The value of the characteristic time f st was then 
calculated for each metal as a function of temperature. 
As an example, the values referred to the liquid tin are 
reported in Table I. 

The values of f st as a function of temperature at 
a total pressure P tot = 10 -1 Pa for a set of selected 


TABLE I Characteristic time t s , as a function of temperature of a liquid tin drop (R = 5 x 10 3 m) from melting to boiling point at a total 
pressure P, M = 10“ 3 Pa and related thermodynamic quantities 


Temperature 

(K) 

Density 
(kg m “ 3 ) 

*0,s 

Pj.s 

(Pa) 

Pj, S 

(Pa) 

Po„, 

(Pa) 

^St 

(s) 

506 

6999.5 

7.20E-10 

j = SnO 

1.20E —18 

j = Sn0 2 
3.90E-25 

4.30E-45 

6.96E + 12 

556 

6966.3 

7.70E-09 

8.30E-16 

3.60E —20 

1.20E-39 

1.14E + 11 

606 

6933.4 

5.60E-08 

1.90E —13 

5.00E-16 

4.20E-35 

3.70E + 09 

656 

6900.8 

3.00E-07 

1.90E-11 

1.60E-12 

2.90E-31 

1.79E + 08 

706 

6868.5 

1.30E-06 

9.90E-10 

1.60E-09 

5.90E-28 

4.28E + 06 

756 

6836.5 

4.50E-06 

3.00E-08 

6.60E-07 

4.30E-25 

4.93E + 04 

806 

6804.8 

1.30E-05 

5.90E-07 

1.30E-04 

1.40E-22 

8.16E + 02 

856 

6773.4 

3.50E-05 

8.20E-06 

1.30E-02 

2.20E-20 

2.16E + 01 

906 

6742.2 

8.40E-05 

8.40E-05 

7.80E-01 

2.00E-18 

8.59E-01 

956 

6711.4 

1.80E-04 

6.80E-04 

3.10E + 01 

1.10E-16 

4.82E-02 

1006 

6680.8 

3.70E-04 

4.40E-03 

8.40E + 02 

4.40E-15 

3.62E-03 

1056 

6650.6 

6.90E-04 

2.40E-02 

1.70E + 04 

1.20E-13 

3.49E-04 

1106 

6620.6 

1.20E-03 

1.10E —01 

2.50E + 05 

2.30E-12 

4.19E-05 

1156 

6590.8 

2.10E-03 

4.40E-01 

3.00E + 06 

3.50E-11 

6.06E-06 


Element: Sn. 

Total pressure: 10“ 3 (Pa). 

Melting point: 505 K. Boiling point: 1158K. 
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TABLE II Values of f sl (s) as a function of temperature at a total pressure P tol = 10 1 Pa for a set of selected liquid metal drops 
(R = 5x 10“ 3 m) 


To. = 10“ 'Pa 

A1 

T m = 933 K 

Ga 

T m = 303 K 

In 

T m = 430K 

Ge 

T m = 1211 K 

Sn 

T m = 505 K 

Pb 

T m = 601 K 

T m 

8.91E + 02 

5.94E + 13 

8.75E + 35 

7.18E—03 

6.96E + 12 

1.52E + 09 

r m +100K 

8.70E + 01 

4.16E + 08 

6.41E + 27 

2.28E-03 

3.70E + 09 

8.21E + 06 

T m + 200 K 

1.32E + 01 

3.70E + 05 

1.90E + 22 

9.15E-04 

4.28E + 06 

1.71E + 03 

T m + 500 K 

- 

1.35E + 01 

8.05E + 12 

- 

3.62E-03 

- 

T m + 600 K 

- 

2.20E + 00 

1.03E + 11 

- 

4.19E-05 

- 



Figure 1 Isotherm trend of the characteristic time t n of liquid tin as a function of the drop radius R at P 10t = 10 “ 3 Pa. (—) characteristic time 

calculated at melting point (505 K); (-) characteristic time calculated at boiling point (1158 K); (—) characteristic time at an intermediate 

T = 870 K; Bold curve: diffusion times calculated for T from 1023 K to 1158K and extrapolated for T from 505 K to 1023 K. 


metals are reported in Table II. By comparing the 
reported values of t sl calculated under fixed conditions 
of temperature (T m < T < T b ) and pressure for differ¬ 
ent liquid metal drops having constant dimensions, 
with the estimated [8-11] corresponding values of 
oxygen diffusion times in the liquid metal (f D = lffi-lO 3 s), 
the two subclasses can be recognized; “small” drops 
when f st »f D and “large” drops when f st «f D , as defined 
by the theoretical model. A tin drop, for example, can 
be considered “small” up to a temperature of about 
700 K, but for temperatures above 1000 K it can be 
considered “large”. Drops of In or Pb can be con¬ 
sidered “small” in the overall range of temperature 
considered here. On the other hand, a drop of Ge can 
be considered always as “large”. 

As a clarifying example, Fig. 1 shows that the classi¬ 
fication in “small” or “large” drops not only depends 
on the drop’s radius but mostly on the liquid metal 
temperature. The trend of f sl as a function of the radius 
of a liquid tin drop is shown at three fixed temper¬ 
atures. The trend of f D is also reported: diffusion times 
have been calculated at the boiling temperature on the 
basis of available literature data [9] and extrapolated 
for T < 1023 K. It can be noticed that despite a small 
radius, at a temperature near the boiling point, 
a liquid tin drop is classified as “large”, whereas for 


a temperature near the melting one the drop is 
“small”. For the intermediate temperatures some 
t st curve may cross the corresponding f D curve and 
a fixed radius can be singled out which discriminate 
the belonging to “small” or “large” drop class. 

7. Conclusions 

A theoretical model to evaluate the interfacial trans¬ 
port of a liquid metal-oxygen system under Knudsen 
conditions has been developed. 

Starting from a local oxygen mass balance in the 
liquid metal drop and taking into account the total 
mass balance, the oxygen molar fraction at the inter¬ 
face and in the liquid phase have been evaluated as 
a function of time. 

Steady-state conditions have been defined for this 
particular system both for the bulk liquid phase and 
the interface. The evaluation of the relevant character¬ 
istic times (for diffusion, and to reach stationary con¬ 
ditions) has allowed two characteristic system “sizes” 
to be defined: “small” and “large” drops, depending on 
the different mechanisms of gas transport in the liquid 
phase. 

The interfacial characteristic times and the trend of 
the oxygen interfacial composition vs. time have been 


311 





calculated for “small” and “large” drops in some 
simple cases. 

Finally, the present model has been used to make 
some estimate of the behaviour of selected metallic 
systems under different conditions. 

An experimental programme is currently under way 
in our laboratory to validate this model; preliminary 
results confirm the reliability of the characteristic 
times to reach a surface steady-state composition esti¬ 
mated by the present work. 
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